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CN : Abstract 

o 

Q . We apply a Gutzwiller-like variational technique to study Joseph- 

"^ ! son conduction across a quantum dot with an odd number of 

electrons connected to two superconducting leads. Our method 
^ ; projects out all states on the dot but the Kondo singlet and 

is valid when Kondo correlations are dominant and no Andreev 
bound states localized at the dot are available for Kondo screen- 
Ch ! ing. In these conditions superconducting pairing is a competing 

fj I effect and the junction is vr— like, to optimize antiferromagnetic 

correlations on the dot. As the superconducting gap increases, 
the Josephson current also increases, but its phase dependence 
becomes strongly non sinusoidal. 
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Phase coherent transport in hybrid superconducting-semiconducting nanos- 
tructures is aheady extensively investigated. In these devices quite interest- 
ing and surprising features emerge, due to electron-electron interaction [1]. 
In this paper we study a Quantum Dot (QD) at Coulomb Blockade (CB) 
with an odd number of electrons N, connected to superconducting contacts 
[2]. Although such a kind of system has not been realized yet, it is very likely 
that continuous improvement of nanofabrication techniques will soon make 
it available. In particular, we address the question whether the formation of 
a Kondo singlet at the dot competes or cooperates with BCS s-wave pairing 
in the leads. We set up a variational non perturbative approach which can be 
adopted as long as A < ksTx, where Tk is the Kondo temperature of the dot 
and A is the superconducting gap. We show how the flow to the Kondo fixed 
point is affected by superconductivity because singlet pairing starts to com- 
pete with Kondo correlations. The system behaves as a vr-Josephson junction 
[3], with a rather small critical Josephson current Ij. As A increases and 
becomes comparable with ksTx, there is a crossover to a regime in which 
Kondo correlations on the dot are very much affected. The junction is still 
vr-like and the Josephson critical current increases but it becomes strongly 
non sinusoidal. 

In our approach bound states localized within the dot area are responsible 
for the bare magnetic moment on the dot. They are only weakly modified 
when the contacts become superconducting, as they turn into Andreev bound 
states within the gap. As discussed in the conclusions, the fact that our 
model does not include other empty Andreev states hinders the possibility 
of screening out the dot spin completely, in the large Kondo limit. 

We believe that this is the reason why the 7r-coupling always wins in our 
model, even when the ratio A/Zc^TV is quite small. Indeed n coupling opti- 
mizes pair tunneling in the presence of strong antiferromagnetic correlations 
on the dot. 

Coulomb Blockade has been widely investigated in Quantum Dots with 
normal leads [4]. If the coupling with the contacts is weak and the charg- 
ing energy is higher than the thermal activation energy, DC conduction is 
strongly dependent on the gate voltage Vg. The DC conductance across the 
dot shows a sequence of peaks, occurring when Vg provides the energy to add 
an extra electron to the dot, taken from the contacts at chemical potential /x. 
As the dot is tuned in a "Coulomb Blockade" valley between two consecutive 
peaks, N is fixed and the conductance is heavily suppressed. 



Kondo conductance may be achieved in a QD with N odd within a CB 
valley by increasing the coupling between dot and leads, provided the tem- 
perature T < Tk [5]. At T < Tft-, a strongly correlated state between dot 
and leads sets in and a resonance in the density of states of the QD opens 
up at n- Correspondingly, DC conductance across the dot increases, until it 
eventually reaches the unitarity limit at T = [6, 7]. 

Since charge dynamics is "frozen out" by CB, the QD in this regime is 
usually modelized as a magnetic impurity with total spin S = ^. A Schrieffer- 
WolfF transformation on the single level Anderson impurity model involving 
virtual zero or double occupancy of the impurity level leads to the effective 
Kondo Hamiltonian with antiferromagnetic coupling between the delocalized 
electrons of the leads and the impurity spin. 

Magnetic impurities embedded in a bulk superconductor are known to 
strongly influence the superconducting critical temperature. Adding impu- 
rity states at energies below the superconducting gap can even give raise to 
gapless superconductivity [8]. 

One can extend all this to the case of superconducting contacts. The 
Schrieffer- Wolff transformation can be performed in the same way, provided 
D » A, where D is the bandwidth of the itinerant electrons. However, the 
resulting model is not equivalent to a single magnetic impurity in a phase 
coherent superconductor. This is a crucial difference with respect to the case 
of normal contacts when the Kondo screening cloud shows that the phase 
coherence is established throughout the system. 

On the contrary, in such a case the contacts retain their individual super- 
conducting properties and their individual phase for the order parameter with 
a phase difference (p. The link between the two superconductors is offered by 
the dot and is tunnel-like. This gives rise to two important features: 

1) Kondo phase coherence and superconducting phase coherence compete; 

2) the system has the properties of a Josephson link between two different 
superconductors. 

The reason is that strong on site Coulomb repulsion makes the phase 
breaking time induced by spin flip scattering processes much shorter than 
the Kondo resonance lifetime h/ksTK [2]. 

Indeed, tunneling across a magnetic impurity between superconductors 
with on site Coulomb repulsion has been recently revisited [9, 3, 10, 11]. If 
A > ksTx, the system is unable to scale toward the strongly-coupled Kondo 
regime. Sub-gap Cooper pair tunneling is strongly suppressed by Coulomb 
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Figure 1: Sketch of cotunneling processes ((a) hole process, (b) particle pro- 
cess) in the perturbative regime, to move a Cooper pair from L to R across 
the QD. The sequence of steps 1 ^ 2 ^ 3 is constrained by Pauli restric- 
tions. Analogous processes not depicted here can take place, in which the 
spins are reversed. 



repulsion, unless each cotunneling step is accompanied by a spin flip at the 
impurity. Pair tunneling can take place but it occurs in a three step sequence 
via virtual empty or doubly occupied intermediate dot states, because of 
Pauli exclusion principle, as depicted in Fig. 1. It has been proposed that 
such a mechanism may reverse the sign of the Josephson current through the 
dot so that the Josephson energy is at a minimum for ip = n, where (f is the 
phase difference between the order parameters of the two superconductors 
attached to the impurity (vr-junction) [3]. This is indeed what we flnd with 
our variational approach. 

Fig. 1 mimicks Cooper pair tunneling via virtual breaking of a pair, one 
at a time, with one quasiparticle moving into and out of the dot. Particle 
and hole processes are allowed across the dot, with a mechanism similar to 
the one involved in the Schrieffer- Wolff transformation. Here we describe 
the hole process only; the particle one is fully analogous. Let us consider 
an I electron localized at the dot. In order to move a pair from L to R, 
we must first empty the dot. The sequence 1 — » 2 — i> 3 corresponds to an 
operator arrangement, CL'\C^p>^CLic'p,p which has a minus sign with respect to 
the standard arrangement for Cooper pair tunneling: c^riCr^CliCli [3]. This 
produces a phase difference of n with respect to the usual Josephson coupling, 
which makes the system a "vr— junction". 



In Superconducting/Normal/Superconducting (SNS) structures, Cooper 
pair subgap tunneling current takes place via Andreev states localized in 
the normal region [12]. When the normal region is given by a QD, charge 
quantization at the dot is not spoiled as the contacts become superconducting 
[13]. The I — V curve of a dot at CB between two superconductors has been 
derived in [14] , showing an interplay between multiple Andreev reflection and 
electron- electron interaction at the QD. 

We expect that flne tuning of the parameters of a S-QD-S device in an 
orthogonal magnetic fleld may allow for the investigation of the full range 
of physical conditions, from A ^ ksTx when the system can flow toward 
the strongly coupled Kondo regime prior to the onset of superconductivity 
in the leads, to A 3> ksTx when perturbation theory holds [11]. In fact, 
the ground state degeneracy required for Kondo coupling to take place can 
be obtained by driving the dot to a level crossing by means of an applied 
magnetic field [7, 15, 16, 17]. 

We study the T = case with A < ksTx using a non perturbative 
variational technique. The fact that a global phase coherence cannot take 
place in the system justifies our variational approach outlined in the following. 
Hence, we add Kondo correlations to a state which is the superposition of 
two BCS states for the left and the right contact with different phases of the 
order parameter. Such a technique has already been applied to the case of a 
dot with normal contacts, and it has been shown to provide good qualitative 
results in the perturbative regime as well, where it reproduces the poor man's 
scaling equation [17]. 

To construct the trial state, we start from the state \(p, s), given by the 
product of the left (L) condensate times the right (i?) condensate, with the 
two order parameters having a phase difference (f, times the state of the dot 
Is): 

|v9,s) = [BCS, L)x|BCS, i?, ^)x|s). (1) 

Here s is the spin component along the quantization axis of the dot spin 
Srf = i at CB with odd A^. 

The minimal model for the Kondo interaction between electrons localized 
on the QD and electrons from the contacts is Hk = J(T{0)-Sd- The spin 
density operator of the delocalized electrons, (t(0), at the position of the dot 



a; = along the vertical axis, is: 
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^(0) = ^E(cl9,a + Ck,q,aK,a'iCL,q',a' + ^R^g'^)- 
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We have used the fermion operators Cgj(Cgj), (j = L,R) in the plane wave 
representation to describe the contacts particles and V is the normalization 
volume of the leads. We take the symmetric case in which hybridization 
of the dot with the L and R contacts, F, is the same. Kondo coupling is 
antiferromagnetic (AF): J > 0. 

The total Hamiltonian is H = Hs + Hk-, where Hs-, defined in eq. (5) 
below, is the Hamiltonian for the L and R superconducting contacts in the 
BCS approximation. 

The correlated trial state is constructed by applying a Gutzwiller-like 
projector Pg to \^,s). As T — > 0, the system scales towards the strongly 
coupled, large J regime. Correspondingly, Pg gradually projects out the 
high-energy components of the trial state, so that eventually only a localized 
spin singlet survives at the QD. 

The "projector" Pg is defined as [17]: 

P,= (l-\g)+g{a{Q)f-^S,-a{Q), 

where g is a, variational parameter which ranges between (7 = and g = 4/3. 
When g = we have Pq = 1 — 48(^-^(0); Pq fully projects out the high energy 
localized spin triplet at x = 0. As g varies from to 4/3 also the localized 
spin doublet state it increasingly projected out. Eventually, when g reaches 
the value g = 4/3, only the localized spin singlet is left over. 
The trial state is defined as: 

\g,v) = Pg\v,s). (2) 

The value oi g{J) is determined by finding the minimum of the energy func- 
tional E[g, J, ip, A], defined as: 

E[g,J,^,A\ = — ■ — . (3) 

{g,v\g,^} 

Eq. (3) can be calculated by first expressing the products PgHPg and (Pg)^ 
in terms of the usual fermion quasiparticle operators aj^q{ajg),(3j^g{(3jg), (j = 

6 



L, R) which destroy (create) excitations on the BCS states of the L and R 
contact and then by normal-ordering the corresponding operator products. 
The operators are [18]: 

Pj,-q = UqC,,-q,i+Vge''^'clg,V (4) 

where Uq and Vg are the BCS coherence factors and (pR = ip while 0^ = 0. 
The Hamiltonian for the contacts Hs is conveniently expressed in terms of 
these operators as: 

Hs = Ebcs+ E Eg{al^a,,q + Pl_gPj,.g). (5) 

q,j=L,R 



Here -Ebcs is the total ground state energy of the condensates and E, 



q 



Q + A^ are the energies of the quasiparticle excitations with ^^ = q^/2m—ix 

(we put h = ks = I throughout the paper). 

The variation of the energy [18] due to the AF coupling at the quantum 
dot in units of the bandwidth D, e[C,,j, 6, ip] = E[g, J, A, fj/D, takes a simple 
form once expressed in terms of the parameters ^ = | (l ^ ffl')) J = 'iJ/^E) 
and 8 = A/D: 



2 (1 + e) + ie - imi + cos(v9))/(2 {N{o)\y 



(l-e + ^^) VTTP-6^cos{^)/N{0)\ 



(1 + e) + ie - imi + cos(v.))/(2 {N{o)\r) ■ ^^^ 

Here A is the BCS electron-electron interaction strength and A^(0) is the 
normal phase density of states at the Fermi level for each spin polarization. 
We will take A^(0)A = 0.3 throughout the paper. The first term is the 
expectation value of Hk, while the second is the raise in the average value 
of the kinetic energy Hamiltonian, Hs, due to the formation of the singlet 
between the QD and the contacts. The value of ^mm, at which the energy is 
at a minimum, measures how much higher-energy spin states are projected 
out: .^inin = corresponds to full projection of states different from a localized 
spin singlet at the impurity. 




Figure 2: e/j vs. C, for different j and 6 = 0.06 in the two cases ip = (broken 
line) and ip = n (full line). For large Kondo coupling j the minimum of the 
curve moves toward ^ = 0. 



At (5 = the strong coupling fixed point is j — > oo and ^min(j) — ^ 0. In 
the flow to the fixed point, there is a large decrease of the Kondo energy: 



fi^ = e[e,j,0,0] 



^n 



(e. 



n-(e. 



(7) 



At (5 7^ 0, we follow the minimum of the energy of the correlated state, 
e*\j,6,ip] = e[^raui,j,S,ip], in the flow j -^ oo. In Fig. 2 we plot e/j vs. 
^ for different j at 6 = 0.06. We see that the minimum of e[C,,j,6,ip] is 
located at smaller and smaller ^ values the larger the Kondo correlations 
are (j -^ oo). However the zero value for ^min, corresponding to the full 
projection of states different from a localized spin singlet on the dot, is never 
reached. Energies corresponding to (f = ii (full line) are always lower than 
those corresponding to (/? = (broken line). This is a general feature of our 
results: the vr-junction behavior is established throughout the whole range 
of 6 values, until A reaches the Kondo scale T^- 

In Fig. 3 we plot ^mm vs. S ioi (p = 0,7i and various j values. In the case 
9? = 0, ^min grows up abruptly for S > 0.25 what shows that the choice ip = 



strongly disfavours the Kondo correlations, as soon as the superconducting 
pairing correlation length ^s decreases. The dependence of ^min on 6 is much 
weaker when (f = tt. Inspection of Fig. 5 shows that two regimes can be 
envisaged as 6 increases: 

a) A -C Tft-: Superconducting pairing in the leads is a perturbation on 
the Kondo correlated state 

For small values of S the minimum of the energy can be attained both for 
ip = and for ip = n. According to Fig. 3, there is a slight decrease of ^min 
when 6 increases for (/? = 0. Hence, the superconducting singlet correlations 
seem to cooperate for small (5's in such a case. By contrast, ^min increases 
steadly with 6 for y? = vr, what shows that singlet correlations on the impurity 
(with characteristic correlation length C,k) compete with singlet pairing but 
are much less disruptive. This is confirmed by Figs. 4a), 4b). The energy 
minimum is strongly affected when 6 increases for (/? = 0, while this is not 
the case for (f = it. Again, the case ip = ir is always favoured in energy for 
any j value (see Fig. 4c)). 

The energy of the correlated state e* [j, 6, ip] , to second order in 6, is: 



^*[^'^'^] - ^-^+^^i^w#-w^°^^' 




2 1 + (e 
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2 iN{o)\y 



+ 

(8) 



Because Q^ < for large j, it is clear that (^ = vr is favoured. 



b) A < Tk- Superconductivity strongly competes with Kondo ordering 
Such a situation is better shown in Fig. 5a): when the phase difference 
between the two superconducting leads is ip = 0, the energy functional, as S 
approaches the limit value 0.3, cannot gain a minimum. Indeed the system 
cannot develop a Kondo singlet on the impurity site. Conversely, Fig. 5b) 
shows that a minimum always exists when the phase difference is ip = it. 
It is easy to derive the Josephson current, according to: 

Ij = 2ede*[j,6,^]/d^. 

In regime a), Ij it is still sinusoidal, of the form Ij = 2e (^(q^x)^ sin ^- 
The sign of such a current is reversed with respect to the conventional one 
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Figure 3: ^min vs. b for various j, for (/? = (broken line) and y? = tt (full line). 
In the case 9? = 0, for small values of 5, ^min is lowered, showing that a small 
amount of superconductivity singlet pairing favors Kondo correlations. On 
the other hand, for large 5, ^ increases abruptly, what signals the disruption of 
the Kondo correlations. By contrast, the weakening of the Kondo correlations 
is much steadier but less sharp when y? = tt. 
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Figure 4: e* vs. j for b = 0.05,0.15,0.29 in the two single cases (/? = (a) 
and (y9 = TT (b) and together (c) for a comparison. 
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Figure 5: e vs. ^ for 5 = 0.05, 0.15, 0.20, 0.22, 0.25, 0.26, 0.27 in the two cases 
(/? = (a) and (/? = vr (b). (a) Deeper minima correspond to smaller 6. For 
6 close to 0.3 the system will not be able to build up the Kondo singlet, (b) 
Deeper minima correspond to smaller 6. 

for Josephson systems, that is, the device behaves as a tt-j unction. In such 
a regime superconductivity starts to compete with Kondo correlations and 
the Josephson critical current comes out to be rather small, as we can see in 
Fig. 6a). 

TT-coupling wins in our model no matter how small the ratio A/ZcbTr- is. 
This is counterintuitive because one expects that, when Kondo correlations 
are fully estabilished and the magnetic moment at the dot is totally screened, 
the scattering of Cooper pairs should be potential-like only. This implies 
that no phase breaking takes place at the tunneling and the junction should 
behave as a conventional Josephson junction, with its energy minimum at 
ip = [10]. This is not the case here because the density of states does not 
include extra subgap Andreev states localized at the dot which could provide 
the screening of the magnetic moment. We expect that, if Andreev bound 
states are accounted for properly in the model, they can provide the required 
full Kondo screening for vanishing A/ksTx ratio. 

As 6 increases, we enter a new regime in which superconductivity strongly 
competes with Kondo correlations. The Josephson critical current becomes 
now sizeable. When 6 exceeds a threshold value (close to 0.3) the Josephson 
current develops strong non sinusoidal components (see Fig. 6b)). 

The regime A > T^, in which the superconducting order is dominant, 
cannot be descibed by our variational Ansatz. In this case, the Kondo inter- 
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Figure 6: Different possible regimes for the Josephson current at zero tem- 
perature for j = 8: a) A ^ Tk] b) A < Tk- 

action is better treated perturbatively. Tliis regime was studied in ref. [10] 
using a NCA perturbative approacli and in ref. [11]. Their resuhs can be 
interpreted in terms of Andreev bound states in the normal region. 

In the noninteracting limit the single particle energy spectrum for the 
S-N-S sandwich has two Andreev levels corresponding to states localized in 
the normal region with energy within the gap. Taking the Fermi level /x as 
the reference energy, the one of the two states with positive (negative) energy 
is mostly particle (hole)-like. In this case, Josephson conduction across the 
normal region involves both Andreev states. Because the levels are non 
degenerate Kondo correlations cannot take place until J is large enough. 
Kondo correlations require that the particle and hole-like states cross each 
other, so that screening of the impurity spin can take place at the dot site. 

In conclusion, we have generalized the variational approach introduced in 
ref. [17] to study Josephson conduction through a quantum dot at Coulomb 
blockade connected to two superconducting leads. Within such an approach, 
we analyzed the region of parameters where the dot lays in the strongly 
coupled Kondo regime and showed that, in this regime, it behaves as a tt- 
j unction with a small Josephson current. Our formalism shows that the onset 
of TT-j unction regime takes place much before antiferromagnetic correlations 
at the dot can be treated perturbatively, that is, much before Kondo effect 
has been disrupted by superconductivity. 
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